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The Asqtad improved staggered fermion formalism has been a valuable tool in successfully cal- 
culating the non-singlet parts of the hadronic spectrum. We are engaged in a project to calculate 
the spectrum of the pseudoscalar singlet mesons with 2 + 1 -flavor Asqtad staggered gauge con- 
figurations. Propagators of flavor-singlet states incorporate contributions from both disconnected 
and connected diagrams, and hence are sensitive to any differences in the actions governing the 
sea and valence fermions on the lattice. As such, they also present the possibility of a probe of 
the validity of the "fourth-root trick" in the staggered fermion formulation. We present an update 
on our progress toward measuring the r\' mass on 2 + 1-flavor Asqtad staggered gauge configura- 
tions, including a review of methods and preliminary results. We also show a strong correlation 
between Tr(/5 <g> 1) and the topological charge in these configurations, as predicted by the index 
theorem. 
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1. Introduction 

The pseudoscalar flavor-singlet meson system is of theoretical interest for a variety of reasons. 
It differs from the non-singlet by the inclusion of disconnected diagrams in its propagator, to which 
is attributed the mass difference between the pion and the tj meson (958 MeV)[jl| §. A thorough 
lattice QCD calculation of the spectrum of the pseudoscalar flavor-singlet system is needed for 
full comprehension of this connection. To date there have been a number of lattice studies of 
Nf = 2 flavor singlet mesons^, ||, ^, [7|, |L [| 10, 11, O]. These proceedings contain a report on an 



N f = 2+1 flavor Wilson lattices fermion study from CPPACS [JT3p. 

We are engaged in a study of the pseudoscalar singlet system with Nf = 2+ I flavors of im- 
proved staggered quarks. The motivation is two-fold. First, the improved staggered "Asqtad" for- 
mulation has a formidable track record of accurate reproduction and prediction of experimentally 
measurable quantities, in part due to the relative high speed of staggered fermion simulation and 
the accessibility of relatively light quark massesjQ]. Second, the staggered formulation has been 
plagued by theoretical questions regarding the validity of the "fourth-root trick" employed in sim- 
ulations where Nf ^ 4. These concerns center on whether or not the det 1 / 4 of the native four-flavor 
staggered fermionic matrix introduces pathologies into the sea quarks. An additional motivation, 
then, is to see whether singlet propagators, which are uniquely sensitive to the fermionic sea, illu- 
minate any such pathologies. 

For Nf degenerate flavors of quarks, the full propagator is 

G n ,(x',x) = (£W)(Y5®l)qitf)Jlqj(x)(Y5®l)qj(x)), (1.1) 

where the 75 1 denotes that the meson has 75 Dirac spinor structure, and is a singlet in staggered 
"taste" space. This expression includes Nf connected terms (pion propagators): 



(L?i(V)(75®l)?/(x')^^(x)(7 5 ®l)^(x)), (1.2) 



and Nf disconnected terms: 



<E?/(^)()S®l)9i(^)E9/W()S®l)9yW>- (1-3) 

i j 

In other words 

G n i(x',x) =NfC(x,x)-NjD(x',x), (1.4) 

where the additional fermion loop in each of the disconnected diagrams gives rise to the negative 
sign. 

As the connected correlators are the propagators for the non-singlet meson we expect, when 
the ground state dominates, for them to decay exponentially in time separation on a Euclidean 
lattice: 

G n (t)=C(t)~e- m « t . (1.5) 

Furthermore, 

G n ,(t) =N f C(t) -NjD(t)^ e - m n' t . (1.6) 
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Then the ratio of disconnected to connected contributions should be 

R(t)= f )' = \-Ae- {m n'- m ^, (1.7) 
V ; NfC(t) 

if the sea quarks are appropriately dynamical. In the case of quenched QCD we expect instead 



[14,15] 



R(t)~A+Bt (1.8) 

It is also possible that the det 1 / 4 in the staggered formulation introduces other pathologies that may 
manifest themselves in this ratio. 

Measuring disconnected correlators is a difficult task. They are noisier than connected corre- 
lators, as they contain fluctuations of the fermionic sea. Hence one must use some sort of all-to-all 
propagators to get as much information out of a particular gauge configuration as possible. Addi- 
tionally, the pseudoscalar singlet is closely related to very slow topological modes, necessitating 
the use of very long time series. We will describe how we address these challenges below. 

2. Simulation and Measurement 

We have at our disposal the large publicly-available library of MILC Asqtad lattices. Our pre- 
liminary study has predominately utilized the "coarse lattice" ensembles, with lattices of 20 3 x 64 



and lattice spacing of a m 0.125fm[|16|]. The typical ensemble has between 400 to 600 configura- 
tions separated by 6 trajectories. 

On each gauge configuration we measure the connected correlators with standard point sources 
with routines integrated into the Chroma code Jl7|]. We measure disconnected correlators with 



stochastic volume-filling sources[18, 19, 20], r\{x). Given that on averaging over sources 

5 v = (T](x) t 7](x)) ?JJ (2.1) 
and defining Mtp = r\ , then the pseudoscalar loop operator is 

*W(0 = L (0 f tonsil? (2-2) 

X4=t 

The operator Ay 5 ®i is the appropriate four-link covariant shift, displacing the quark and antiquark 
to opposite corners of the hypercube, with appropriate Kogut-Susskind phasing to effect the 75 (g> 1 
operator. We average over N s = 64 noise sources. The disconnected correlator is then: 

£>(0 = 7l>('W+0 (2-3) 
1 f 



In practice one can do better than Equation 2.2. Venkataraman and Kilcup[10] describe a 
variance reduction (VKVR) trick applicable to the four-link pseudoscalar singlet operator. Recog- 
nizing that the Asqtad Tp connects only odd sites to even sites, and the M^M = (D 2 + m 2 ) and its 
inverse connects even to even and odd to odd only then: 

mi^Ays&Q) = (^A K0 iT)), (2.4) 
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with the left-hand expression having a minimized variance. The equivalence of the expectation 
value depends on the A^gi operator separating the source and sink by an even number of links 
(four). 

The numerical work by the TrinLat group suggests that using a set of dilute noise sources, 
each defined on a subset of the lattice but collectively spanning the lattice, can reduce the variance 
of measured propagators^]. In the case of disconnected correlators for pseudoscalar singlets with 
Asqtad fermions, we found no discernible advantage. 



□ Nf=2+lbeta=6.76am=0.010.0.05 658 contigs 
Nf=2+lbeta-6.76am=0.007,0.05 437 configs 
O Nf=0 beta=8.00 am=0.05 4819 configs 



Figure 1: A comparison of D/C ratios for quenched and Nf = 2 + 1 flavor configurations. The /3 = 2.76 
am = 0.007 points were computed with only 40 noise sources per configuration and without the VKVR trick. 









— MILC ensemble: 408 cfgs 

— stream 0: 408 configs 
stream 1 : 408 configs 

stream 2: 408 configs 

— stream 3: 408 configs 
stream 4: 408 configs 

stream 5: 408 configs 

stream 6: 408 configs 

stream 7: 408 conl'i«s 

stream 8: 408 configs 

stream 9: 408 conl'i«s 

all 4819 configs 
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Figure 2: D/C ratio for 4819 j3 = 8.00, am = 0.05, quenched configurations, and for subsets of 408 
configurations, including the original MILC ensemble. 



3. Analysis 

We were able to extract signals for connected and disconnected correlators for several ensem- 
bles of MILC coarse lattices. For Nf = 2 + 1 configurations with two degenerate flavors of light 
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Figure 3: Error on R(t = 8) for 4819 /3 = 8.00, am = 0.05, quenched configurations, as a function of 
jackknife bin size. 



sea quarks and one strange sea quark flavor, expression L7 generalizes to[pi|] 



AD qq {t)+AD qs {t)+D ss {t) 

R{t) ~ 2c M (0+c M (0 • (3 " 1} 

The disconnected correlators, D qq , D qs and D ss , are formed with two light quark loop operators, 
with a light loop operator and a strange loop operator, and with two strange loop operators, re- 
spectively. The connected correlators C qq and C ss are light quark and strange quark correlators, 
respectively. 

When we form the D/C ratio for Nf = 2+1 lattices we get a curve marginally consistent with 
a plateau at one as in Equation |L7[ However, it is also consistent with a linear relationship such 



as that in Equation 1.8. With the MILC quenched ensemble (j3 = 8.0 am=0.050, 408 configs), we 



find a D/C curve that is inconsistent with a linear relationship as in Equation L8 See Figure [I[ 



We hypothesized that the 75 (g> 1 loop operator might be subject to the particularly long au- 
tocorrelation times of slow topological charge modes. As quenched configurations are relatively 
inexpensive to produce, we have extended the MILC ensemble by nearly a factor of 12 to 4819 con- 
figurations. When we analyze the entire ensemble, the R(t) curve is more consistent with Equation 



L8| (see figure |2|). However, subsets of 408 configurations are often mutually inconsistent by a 
margin of two standard deviations even when determined with relatively large jackknife bin size 
— far greater than the integrated autocorrelation time of about 10 configurations. We interpret this 
to mean that the time series autocorrelation does not decay as a single exponential, and that there 
are autocorrelation time scales in this mode that are as long as 100 configurations or more. 



4. Topological Charge 



Using the Atiyah-Singer index theorem [22], Smit and Vink [ 23 ] equate a fermionic expression 
for the topological charge 



Q = mK p {Tr{r 5 M- i )) l 



(4.1) 
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with the traditional gluonic definition: 

Q(x) = 



64k 2 



(4.2) 



Equation ^Jj * s merely the loop operator in Equation [Z^, integrated over time. In [|24[], Alles et at 
show the equivalence of these two definitions on Wilson configurations. As a cross check of our 
loop operator, we measured the topological charge using the gluonic definition, Equation A_2 , using 
hypercubic blocking^] implemented in the MILC code [26]. We found a strong linear correlation 
between the integrated loop operator and the topological charge as defined gluonically in Equation 
4.2. (See Figure EJ for an example.) 




Figure 4: Fermionic versus gluonic topological charge (after three HYP cooling sweeps) for Nf = 2+1 
j3 =6.76 am = 0.01 , 0.05. Normalizations not consistent between the axes. 



5. Conclusions 

We interpret these results to mean that there is potential in Asqtad fermions as an avenue 
for investigating pseudoscalar singlet systems. If our experience with quenched configurations is a 
guide, one can achieve good signals for disconnected propagators with the order of 10 4 trajectories. 
While this is far in excess of any current MILC Asqtad ensembles, we are currently in the process 
of generating such a ~ 2 x 10 4 trajectory ensemble on two racks of a QCDOC machine. Use of the 
VKVR trick has proven to be invaluable for reducing the variance of disconnected 75 1 operators. 
Careful measurement of a suite of fuzzed operators may also help in extracting the pseudoscalar 
singlet system spectrum for Nf = 2 + 1 staggered QCD, as well as shed light on the validity of the 
fourth-root trick. 
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